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Sufficiently inclusive observables in the decay of the tau lepton can be calculated using the methods of per- 
turbative QCD. These include the asymmetry parameter A T that determines that angular distribution of the 
total hadron momentum in the decay of a polarized tau. It should be possible to measure A T accurately using 
existing data from LEP. Reliable estimates of theoretical errors are essential in order to determine whether a given 
observable is sufficiently inclusive to be calculated using perturbative methods. The theoretical uncertainties due 
to higher orders in a 3 can be estimated using recent calculations to all orders in the large-(33 — 2Nf ) limit. These 
estimates indicate that tau decay data can be used to determine a s (Mz) to a precision of 2% or better. 



1. INTRODUCTION 

What can perturbative QCD tell us about the 
decays of the tau lepton? If you ask this ques- 
tion to the typical man on the street, he will an- 
swer "Absolutely nothing!" Perturbative QCD 
tells us about the interactions of quarks and glu- 
ons with large momentum transfer. But the decay 
of the r is dominated by decays into single par- 
ticles and resonances: 7r, p, a%, etc. The QCD 
interactions that bind quarks and gluons into 
these hadrons necessarily involve small momen- 
tum transfers and are therefore completely out- 
side the domain of perturbation theory. There- 
fore one should not expect perturbative QCD to 
tell us anything about the decays of the r. 

But suppose we ignore this objection, and pro- 
ceed to calculate the decay rate of the tau into 
a neutrino plus quarks and gluons using per- 
turbative QCD. The resulting expression for the 
hadronic decay rate is an expansion in powers of 
a s . Using the measured value of the hadronic 
decay rate, we can determine a s . We find that 
this value for the QCD coupling constant agrees 
amazingly well with the best determinations of 
a s . Is this just a remarkable coincidence? Or is 
it possible that sufficiently inclusive observables 
in tau decay really can be calculated using the 
methods of perturbative QCD? 



* Invited talk presented at the TAU 96 Workshop, Estcs 
Park, Colorado, September 1996. 



2. QCD PREDICTION FOR R T 

The QCD prediction for the ratio R T of the 
hadronic and electronic branching fractions of the 
tau has the form [Dj 



Rt = 3 {\V ud \ 2 + \V us \ 2 ) S EW 
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The electroweak corrections include a multiplica- 
tive factor 5ew = 1-0194 and a small additive 
correction Sew — 0.0010 0]. There is a perturba- 
tive QCD corrrection S per t that can be expressed 
as an expansion in powers of a s at the scale M T 
(or any other scale /i of your choosing) H: 
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The remaining QCD corrections are the power 
corrections 5 power , the most important of which 
are the following B: 



-16 



A/2 
ll^ 2 



32tH 



2 



«.(M r )V j^GG) 

(ipTipipTip) 

AC? ' 



terms of the form 



(3) 



where m 2 is a weighted average of the running 
quark masses m 2 , m 2 , and m 2 evaluated at 
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the scale M T with weights 1/2, |T4d| 2 /2, and 
|K, s | 2 /2, respectively. Similarly, (mtpijj) is the 
same weighted average of the quark condensate 
matrix elements (m u uu), {vriddd), and (m s ss). 
The matrix element (^-GG) is called the gluon 
condensate. 

The QCD prediction (1) is based on the fact 
that the inclusive hadronic decay rate of the tau 
involves only 2 momentum scales: the mass M T 
and Kqcd, the scale associated with nonpertur- 
bative effects in QCD. One can imagine increas- 
ing M T while holding the QCD coupling constant, 
and therefore Aqcd, fixed, ft is possible to sys- 
tematically separate the effects of "hard" partons, 
whose momenta scale with M T , from the effects 
of "soft" partons, whose momenta remain propor- 
tional to Aqcd as M T increases. At very short 
distances of order 1/M\y, the decay of the t~ 
proceeds through the decay into a neutrino plus 
du (or su). The QCD corrections involve both 
hard partons and soft partons. The perturbative 
correction S per t consists of corrections from the 
emission of hard partons and from the exchange 
of hard partons. It can be calculated as a pertur- 
bation series in a s (M T ). The first 3 terms in this 
series are known and they are given in (2). 

What about the effects of soft partons? There 
are large corrections from the emission of soft 
gluons and large corrections from the exchange 
of soft gluons, but they cancel order-by-order in 
a s . In other words, the dominant effects of the 
soft gluons can be expressed as a unitary trans- 
formation that does not change the decay rate. 
At the perturbative level, this unitary transfor- 
mation describes the evolution of an initial state 
consisting of hard partons emerging from the de- 
cay into final states consisting of hard partons 
and soft partons. At the nonperturbative level, 
the soft gluons have dramatic effects, binding the 
hard partons into final state hadrons. Neverthe- 
less, the dominant effect of the soft gluons can still 
be expressed as a unitary transformation which 
evolves an initial state consisting of hard partons 
into final states consisting of hadrons. Since a 
unitary transformation preserves probabilities, it 
does not change the decay rate. 

Now the effects of the soft partons are not ex- 
actly a unitary transformation. There are cor- 



rections which can be taken into account through 
power corrections S power , like those given in (3). 
The term involving the gluon condensate (-GG) 
takes into account effects of soft gluons that can 
not be expressed as a unitary transformation. 
The quark condensate (miptl>) takes into account 
corrections from soft quarks and antiquarks that 
are suppressed by the quark mass. The matrix el- 
ements of the form ('t/jTi/j^Ti/j) take into account 
corrections from soft quarks and antiquarks that 
would be present even if the quarks were massless. 

The QCD predictions for all inclusive observ- 
ables in tau decay have the same general form 
as the predictions for R T in (1). There is a free- 
quark prediction that corresponds to the decay 
into a neutrino plus du or su. There are per- 
turbative corrections that can be expressed as an 
expansion in a s (M T ) like that in (2). Finally, 
there are power corrections like those given in 
(3). They include corrections from running quark 
masses and soft parton corrections that are ex- 
pressed in terms of vacuum matrix elements. 

3. OBSERVABLES 

Below, I enumerate the observables in tau de- 
cay that can be calculated using the methods of 
perturbative QCD. 

3.1. Inclusive decay rate 

The QCD prediction for the inclusive hadronic 
decay rate normalized to the electronic decay is 
given in (I). The ratio R T can be further resolved 
into the contributions from the weak vector cur- 
rent and the weak axial- vector current. It can also 
be resolved into "non-strange" contributions pro- 
portional to | V^d | 2 and "strange" contributions 
proportional to \V US \ 2 

3.2. Invariant-mass distribution 

Let s be the square of the invariant mass of the 
hadrons in the decay of the r. The distribution 
dR T /ds, or more accurately the moments of this 
distribution, can be calculated using perturbative 
methods. The th moment is simply R T itself. 
The QCD predictions for the moments are given 
in Ref. ||. These moments have been used to 
measure the matrix elements that appear in the 
power corrections (3) ||. 



3 



3.3. Angular distribution 

There is one inclusive observable in tau decay 
that has not yet been measured experimentally. 
This is a parameter A T that determines the angu- 
lar distribution of the total hadron momentum in 
the decay of a polarized tau. Let 9 be the angle in 
the tau rest frame between the total hadron mo- 
mentum and the quantization axis for the spin 
of the r. If the tau is unpolarized, the angular 
distribution is uniform in cos 9. If the tau has 
polarization P, the angular distribution is 

^^ = ip T (l+A T Pcos0). (4) 
a cos 9 2 

The size of the asymmetry in cos 9 is determined 
by the parameter A T . This parameter is known 
for various exclusive decay modes. It has the 
value +1 if the hadronic final state is a single 
7T. The value is close to —1 for the transverse po- 
larization modes of the p and a\, and close to +1 
for the longitudinal polarization mode. The in- 
clusive asymmetry parameter A T is the weighted 
average over all hadronic final states of the asym- 
metry parameter for each of the exclusive decay 
modes. 

The asymmetry parameter A T can be calcu- 
lated using the mothods of perturbative QCD 
with the same degree of rigor as R T itself . The 
free-quark prediction for A T is |. The QCD cor- 
rections increase the prediction to A T = 0.415 ± 
0.022. Thus the contributions of the n, p, a\ and 
all the multihadronic modes must conspire to give 
a value that is about 25% larger than the free 
quark value. 

An "experimental" value for A T has recently 
been obtained by Weinstein (|]. He generated 
50,000 polarized tau decay events by running the 
Monte Carlo programs KORALB and TAUOLA 
and obtained the value A T = 0.385±0.007. These 
programs use the known branching fractions for 
the various hadronic decay modes. They incorpo- 
rate the correct asymmetry parameters for the 7r, 
p, and a\ modes, but they are not tuned to give 
the correct asymmetries for the multihadronic 
modes. Thus a real experimental measurement 
of A T is needed in order to test the QCD predic- 
tion. 



A good way to measure A T is to use the taus 
from Z° decay, because they are naturally polar- 
ized with polarization P = —0.14. Suppose one 
could assemble an unbiased sample of taus from 
Z° decay. If, for each tau decay, one measured 
the invariant mass s of the hadrons and their to- 
tal energy E in the Z° rest frame, then A T could 
be determined by calculating the following aver- 
age over the r sample: 

From simple statistics, one would need a sam- 
ple of at least 2000 taus to distinguish the free 
quark prediction 0.508 for (5), which corresponds 
to A T — i, from the unpolarized value of 0.5. 
One would need at least 30,000 taus to distinguish 
the QCD prediction 0.5097± 0.0005 from the free 
quark prediction. With a sample of 400,000 taus 
from Z° decay, it should be possible to measure 
A T rather accurately. Note that the prediction 
that the asymmetry parameter for all the exclu- 
sive modes must average out to a value near the 
free quark value i is already rather remarkable 
and thus even a crude measurement of A T would 
be useful. It would be the first confirmation of a 
perturbative QCD prediction for spin-dependent 
observables in tau decay. An accurate measure- 
ment of A T in agreement with the QCD predic- 
tion would provide dramatic evidence that inclu- 
sive observables can be accurately calculated us- 
ing the methods of perturbative QCD. It should 
lay to rest any questions of whether the accuracy 
of the value of a s obtained from r decay is merely 
fortuitous. 

4. ERROR ESTIMATES 

In order to make quantitative tests of the QCD 
predictions for the observables described in the 
last section, it is essential to have reliable esti- 
mates of the theoretical errors. In the proceed- 
ings of TAU 94, both Altarelli and Narison made 
attempts to estimate the theoretical errors in the 
determination of a s (Mz) from r decay. Their es- 
timates are given in Table 1. The bottom line is 
that Altarelli's estimate of the theoretical error is 
larger than Narison's by a factor of 3. This fac- 
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running: M T — > M z 


20 


10 


freezing of a s 


10 




quark masses 




5 


condensates 




9 


higher orders in a s 


65 


14 


fx dependence 




9 


renormalization scheme 




5 


"other theoretical errors" 


10 




Total Error 


70 


23 



tor of 3 makes all the difference in the world. If 
Narison's estimate is correct, then tau decay re- 
mains competitive with the best determinations 
of a s that are available. If Altarelli's estimate is 
correct, then the determination of a s from tau 
decay is only interesting as a qualitative test of 
perturbative QCD. 

Before assessing the error estimates of Altarelli 
and Narison, it is necessary to discuss some of 
the theoretical issues that are involved in these 
estimates. 

4.1. Renormalization scale /x 

The ratio R T , if calculated to all orders in 
a s , would be independent of the renormalization 
scale fi. However, if the perturbation series is 
truncated at some order in a s , it will depend on 
fi. The value of a s determined by measuring R T 
will therefore depend on the choice for /j. 

This can be illustrated by considering the per- 
turbation series truncated after the order— a s 
term. If we choose the renormalization scale 
to be (i = M T , the prediction for R T (ignor- 
ing electroweak and power corrections and setting 



\V ud \ 2 + IK.. 



3 H 



2 = 1 for simplicity) is 
a s (M T 



(6) 



If we choose the scale to be some fraction x of the 
tau mass \i = xM r , then the prediction is 

a s (xM T 



R 



3 1 + 



(7) 



Having determined a s (xM T ) from the measured 



value of R T , we can then determine a s (M T ) by 
using the renormalization group: 



a s (M T ) 



a s (xM T ) 



1 + £a.(xM T )log(l/a;)' 



(8) 



The resulting numerical value of a s (M T ) will dif- 
fer from that obtained directly from (6). The 
difference between (6) and the combination of (7) 
and (8) amounts to summing certain terms of the 
form a™ log n (M T //z) to all orders in n. 

The dependence on the renormalization scale fj, 
decreases if the perturbation series is calculated 
to higher order in a s . However, as long as the 
perturbation series is truncated, there will always 
be some dependence on fi. 

4.2. Truncation of the perturbation series 

In (2), we have truncated the perturbation ex- 
pansion for R T in powers of a s (M T ), but this is 
not the only way to truncate the perturbation se- 
ries. The ratio R T can be expressed in the form 
of a contour integral 

.„ o i r dt 



12tt z 
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1 - J H? + a BrB?) I '("-<- t >), (9) 



where the function D(a s ) has a perturbation ex- 
pansion in powers of a s : 



D{a s ) 
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The contour integral in (9) runs counterclock- 
wise around the circle |i| = M 2 , beginning at 
t = Af 2 + ie and ending at t = M 2 - ie. The 
function a s (— t) in (9) is the analytic continu- 
ation of the running coupling constant to the 
complex i-plane. It is defined by the bound- 
ary condition a s {— t) = a s (M T ) at t = — M 2 
and by the differential equation t(d/dt)a s (—t) = 
f3(a s (—t)), where f3(a s ) is the beta function for 
QCD: /3(a s ) = -(9/47r)af + .... 

One way to truncate the perturbation series is 
to first carry out the contour integral and then 
truncate the expansion for R T . This gives the 
usual perturbation series for R T given in (2). An 
alternative possibility is to truncate the expan- 
sion for D(a s ) and then carry out the contour 
integral analytically. This has been advocated in 
particular by deLiberder and Pich ||. The re- 
sulting expression for R T is 



Rr = 3 1 + 04± + iJ *™ 
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where the coupling constant with the bar over it 
represents the following weighted average of its 
values along the contour in the complex plane: 

1 fin 

a n s {-t) = 7^ d9 (1 - 2e 10 + 2e m - e ue ) 

xa"(-M r V e ). (12) 

The effect of this prescription for calculating 
R T is more easily illustrated using the simpler 
example of the ratio R e + e -. It can also be ex- 
pressed as a contour integral involving the same 
function D{a s ) as in (9): 

iWU/iHl27r 2 -^ / ±D(a s (-t)). (13) 
2tt« J\ t \=s t 

We choose the center-of-mass energy a/s of the 
e + e~ to be below charm threshold, so that the 
number of light quark flavors is 3, just as in tau 
decay. If we integrate and then truncate R e + e -, 
then at one loop, we obtain 



iWU/i) = 3 {i + ^} 



(14) 



If we truncate D(a s ) and then integrate, we ob- 
tain 



R 



!+ e-(Vi)=3 1 + ^ , 



(15) 



where a s (—t) is a s (—se ) averaged over the an- 
gles 8. Expressed in terms of a s (y / i), this average 
is | arctan (|a s (y / s)) . Its expansion in powers of 



a s (-t) 



1 



27tH 
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Thus the expression (15) for R e + e - differs from 
(14) by the resummation of terms of the form 
7r 2 "(a s /7r) 2n to all orders in n. 

Similarly, the perturbative expansion (11) for 
R T differs from the one given in (2) by the re- 
summation of terms of the form 7r 2 ™(a s /7r) 2 " to 
all orders in n. DeLiberder and Pich give a num- 
ber of arguments why this resummation should 
be preferred, but I don't find them convincing. 
One argument is that it greatly decreases the de- 
pendence of R T on the renormalization scale \x. 
This is true, but there is an equally convincing 
argument against the resummation, and that is 
that the power series expansion for D(a s ) is more 
severely divergent than that for R(a s ). 

4.3. Divergence of the perturbation series 

It has been known for a long time that most 
perturbation expansions in field theory, includ- 
ing QED and QCD, are actually divergent series. 
Suppose you were able to calculate the correction 
of order a™ for any value of n. For small val- 
ues of n, you might find that adding more and 
more correction terms gives a better and bet- 
ter approximation. But if you continued adding 
higher and higher orders in a s , you would eventu- 
ally find that the approximation would get worse 
and worse. The order in a s at which the series 
begins to diverge depends on the value of a s , de- 
creasing roughly as l/a s - Thus the divergence 
of the perturbation series is a much more impor- 
tant issue for tau decay than it is for applications 
of perturbative QCD at higher energy, where the 
running couping constant is smaller. 
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There has been a dramatic development in the 
last five years that has made the problem of the 
divergence of the perturbation series much more 
concrete. A subset of the Feynman diagrams for 
R T and i? e + e - have been calculated to all orders 
in a s |12|| . The diagrams are those of order a™ 
with the maximum number n — 1 of quark loops. 
The complete perturbation series for R T has the 
form 

Rr = 3 jl + + f> n <(M T ) j . (17) 

The coefficients r n are polynomials in the number 
Nf of light quark flavors, or equivalently, in 33 — 
2N f : 



+ r^(33-2JV» + .. 
+ ri"- 1 >(33-2iV / ) r 



(18) 



It is the numbers r„ that have been computed 
to all orders in n. 

Given the coefficients ri™ -1 -' defined in (18), we 
can calculate a subset of the higher order correc- 
tions to R T to all orders in a s . I will denote the 
resulting expression by R T and refer to it as the 
large - (33 - 2N f ) limit: 

K = 3 jl + + pj n a«(M T )\ , (19) 

where f„ = rl™ _1) (33 - 2N f ) n - 1 . The first few 
terms in the expansion for R T are 



R T 



3 1 



28i 
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a s (M T 



a s (M T ) 



a s (M T ) 



5.1 



a s {M T ) 



156.7 



a s (M T ) 



(20) 



This expansion can be compared with the exact 
expansion for the corrections to R T , which is given 
in (2). Note the remarkable agreement between 
the coefficients of a 2 s and a 3 s in the two expan- 
sions. This gives us some confidence that the 
large-(33 — 2Nf) limit can predict correctly the 
sign and order of magnitude of the coefficient of 
the at term. 



Since we know the expansion (19) for R T in the 
large~(33 — 2Nf) limit to all orders in a s , we can 
use it to study the behavior of the perturbation 
expansion |l3|| . One finds that (19) is a divergent 
series with coefficients f n that grow asymptoti- 
cally like n\: 



-53 
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Thus the sum in (19) if taken literally is meaning- 
less. However, there is a standard procedure for 
recovering an analytic function from the divergent 
power series generated by its Taylor expansion. It 
is called Borel resummation. From the series (19) 
for R T , one constructs the Borel transform B(b) 
by dividing each coefficient by (n— 1)!: 



B(b) = 



n— 1 v ' 



(22) 



The power series for B(b) is less divergent than 
that for R T , and actually converges for \b\ < 4ir/9. 

An analytic expression for B(b) is known. It 
has poles at integer multiples of 4ir/9. There are 
poles on the negative real axis at b = n(4%/9), 
n = — 1, — 2, . . . that are called ultraviolet renor- 
malons. There are also poles on the positive real 
axis at b = n(47r/9), n = 2, 3, . . . that are called 
infrared renormalons. Given the analytic expres- 
sion for the function B(b) defined by (22), we can 
recover the desired function R T by computing the 
inverse Borel transform: 



Rr = 3 i 1 



db exp 



a s (M T ) 



B(b) . (23) 



The infrared renormalon poles on the integration 
contour can be handled using a principal value 
prescription. 

Now the divergence of the perturbation series 
for R T is dominated by the singularity of B(b) 
that is closest to the origin. This is the first ultra- 
violet renormalon, which is a pole at b = —4tt/9. 
The behavior of B{b) near this pole is 



B(b) — > —t 

157T 



-5/3. 



1 



(24) 



4tt 



The remainder is analytic inside the circle |6| < 
87r/9. The pole term can be expanded as a power 
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series in b, but the series has a radius of conver- 
gence of only 4n/9. Thus perturbation theory 
gives a convergent approximation to the function 
B{b) in the integrand of (23) only in the interval 
< b < Att/9. Outside that interval, the pole 
term given in (24) cannot be approximated by a 
truncated perturbation series. The integral of the 
pole term from 47r/9 to infinity therefore gives a 
"renormalon correction" to R T /3: 

I db exp 

Utt/9 

1 

X - 



-5/3 



15tt 



a s (M T ) 



1 + 



(25) 



The asymptotic value of the integral for small val- 
ues of a s {M T ) is 



1 

15« 



-5/3 a s( M r) 



exp - 



47T 



(26) 



9a s (M T ) / 

Finally, using the expression 27r/(91og (M T /A)) 
for the running coupling constant a s (M T ), the 
leading renormalon correction reduces to 



± e - 5 /3^s{M T ) A 2 



15 - tt Mr (27) 

where A is the renormalization group invariant 
scale parameter in the MS scheme. 

The analysis above was carried out for the 
standard truncation of the perturbation series 
for R T . With dcLibcrder-Pich rcsummation, one 
must apply Borel resummation to the function 
D(a s ) in the integrand of (9). The Borel trans- 
form of this function has a double pole at b = 
— 47r/9. The renormalon correction analogous to 
(25) therefore has two terms corresponding to the 
double pole and a single pole. After integrating 
over t, it reduces to 

dbcxp ( 

/4tt/9 V 

2 5 

+ 



1. e -5/3 [ C 
9?r J 47l 



a s (M T ) 



(1 + ^) 2 



1 + 



12sin(§6) 



f*(l-&6)(3-&&)(4-&*) 



(28) 



The asymptotic value of the integral for small val- 
ues of a s (M T ) is 

_ 4 5/3 a s (M T ) A 2 



Mf 



(29) 



Note that this is larger by a factor of 20 than 
the corresponding correction (27) for the stan- 
dard truncation. This simply reflects the fact that 
the perturbation scries for D(a s ) diverges more 
severely than that for R T . 

4.4. Estimates of the perturbative error 

We are now in a position to assess the error 
estimates of Altarelli and Narison. The pertur- 
bation series for R T can be written 



R T = 3|l + 5^r„a?(M T ) + r 4 a 4 (M T 

^ n=l 

OO N 

+ £r n a?(M T ) . 

n=5 



(30) 



The first 3 correction terms in the expansion are 
given in (2). We have separated the unknown a 4 
term from the sum of all the higher-order correc- 
tions. 

Narison's error estimate is based on the as- 
sumption that the error from the a 4 term in (30) 
dominates over that from the sum of all higher 
orders. He used the values of the lower-order co- 
efficients n , T2 , and r3 to guess a reasonable range 
for the values of r 4 . The resulting estimate of the 
error from higher orders is 



(Apert)Na 



±50 



a s (M T ) 



(31) 



Translating this into an error on a s (Mz), we ob- 
tain the value ±14 x 10~ 4 given in Table 1. 

Altarelli's estimate was based on the assump- 
tion that the error from the sum of all higher 
order terms in (30) dominates over that from the 
a\ term. As shown in the previous subsection, 
there is a renormalon contribution to that sum 
that is proportional to A 2 /M 2 . Altarelli's error 
estimate is 



(Ap e rt) Altarelli — ± 



1 A2_ 

45?' 



(32) 



Translating this into an error on a s (Mz), we ob- 
tain the value ±65 x 10~ 4 given in Table 1. There 
is no apparent calculation underlying the coeffi- 
cient 1/4 in (32). It seems to be pure guesswork 
on the part of Altarelli. 
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Table 2 

Estimates of the theoretical error in a s (M T ) (in units of 10 -4 ) from higher orders in a s (M T ). 

error from error from sum 

n = 4 term of n > 5 terms 



large-(33 - 2N f ) limit (truncation in a s {M T )) 40 0.2 

large-(33 — 27V/) limit (deLiberder-Pich rcsummation) 4 4 

Altarelli 65 

Narison 14 



One of the problems with both of these error 
estimates is that they require guessing the mag- 
nitude of unknown coefficients: r± in the case of 
Narison, the coefficient of A 2 /M^ from the renor- 
malon correction in the case of Altarelli. One way 
to avoid such guesswork is to use the explicit cal- 
culations to all orders in the large-(33 — 27V/) to 
estimate the errors. As an estimate of the error 
from the a\ term in (30), we can use the magni- 
tude of the a\ term in the large-(33 — 27V/) limit: 



A pert = ±\f 4 \a 4 s {M T ). 



(33) 



As an estimate of the error from the sum of all 
higher order terms in (30), we can use the mag- 
nitude of the sum of all higher order terms in the 
large-(33 - 27V/) limit: 



^pert 



± 



^f„a"(M T 



(34) 



The divergent series on the right side of (34) is 
defined by Borel resummation. A simpler esti- 
mate, which can be compared directly with that 
of Altarelli, is to pick out the leading renormalon 
correction, which is given in (27) for the standard 
truncation in a s (M T ) and in (29) for deLiberder- 
Pich resummation. 

The error estimates based on the large-(33 — 
27V/) limit are given in Table 2. The standard 
truncation gives a much larger error from the a 4 
term and a much smaller error from the sum of all 
higher-order terms. With deLiberder-Pich resum- 
mation, there is a better balance between these 
two sources of uncertainty. Adding the two errors, 
we take the total error in a s (Mz) from higher or- 
ders in cn s to be 8 x 10~ 4 . The error estimates of 



Altarelli and Narison are also given in Table 2 for 
comparison. The error estimate from the large- 
(33— 2Nf) limit is consistent with that of Narison, 
but much smaller than that of Altarelli. The ex- 
plicit calculations in the large-(33 — 27V/) limit 
indicate that the coefficient of the A 2 /M 2 correc- 
tion arising from the first ultraviolet renormalon 
is an order of magnitude smaller than assumed 
by Altarelli. 



5. CONCLUSIONS 

Sufficiently inclusive observables in tau decay 
can be calculated using the methods of perturba- 
tive QCD. A particularly interesting observable 
that has yet to be measured is the asymmetry 
parameter A T that describes the angular distri- 
bution of the total hadron momentum in the de- 
cay of a polarized tau. Using existing data from 
LEP, it should be possible to measure this param- 
eter with sufficient precision to discriminate be- 
tween the free-quark prediction of 3 and the QCD 
prediction A T = 0.41 ± 0.02. Such a measure- 
ment would be the first test of QCD predictions 
for spin-dependent observables in tau decay and 
it would demonstrate conclusively that methods 
based on perturbative QCD give accurate predic- 
tions for sufficiently inclusive observables. 

In order to determine whether a given observ- 
able in tau decay is sufficiently inclusive to cal- 
culate using perturbative methods, it is essential 
to have reliable estimates of the theoretical er- 
rors. One can use recent calculations to all or- 
ders in a s in the large-(33 — 27V/) limit to esti- 
mate the errors from the higher-order terms in 
the perturbation series. The resulting error esti- 



mate on a s (M z ) is consistent with the 1994 esti- Nucl. Phys. B463 (1996) 511. 

mate of Narison, but much smaller than the error 
estimate of Altarelli. These explicit calculations 
indicate that Altarelli overestimated by an order 
of magnitude the size of A 2 /M 2 corrections asso- 
ciated with the first ultraviolet renormalon. The 
bottom line is that a s {Mz) can indeed be deter- 
mined from T-decay data with a precision of 2% 
(or perhaps even better). 

This work was supported in part by the U.S. 
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Physics, under Grant DE-FG02-91-ER40684. 

REFERENCES 

1. E. Braaten, S. Narison, and A. Pich, Nucl. 
Phys. B373 (1992) 581. 

2. W.J. Marciano and A. Sirlin, Phys. Rev. Lett. 
61 (1988) 1815; E. Braaten and C.S. Li, Phys. 
Rev. D42 (1990) 3988. 

3. K. Schilcher and M.D. Tran, Phys. Rev. D29 
(1970) 570; E. Braaten, Phys. Rev. D39 
(1989) 1458. 

4. E. Braaten, Phys. Rev. Lett 60 (1988) 1606; 
S. Narison and A. Pich, Phys. Lett. B211 
(1988) 183. 

5. F. deLiberder and A. Pich, Phys. Lett. B289 
(1992) 165. 

6. L. Duflot, Nuclear Phys. B (Proc. Suppl.) 40 
(1995) 37; A. Hoecker, proceedings of TAU 
96, Estes Park, Colorado (Sept. 1996). 

7. E. Braaten, Phys. Rev. Lett. 71 (1993) 1316. 

8. A. Weinstein, private communication. 

9. G. Altarelli, Nucl. Phys. B (Proc. Suppl.) 40 
(1995) 59. 

10. S. Narison, Nucl. Phys. B (Proc. Suppl.) 40 
(1995) 47. 

11. F. deLiberder and A. Pich, Phys. Lett. B286 
(1992) 147. 

12. M. Beneke, Nucl. Phys. B405 (1993) 424; 
D.J. Broadhurst, Z. Phys. C58 (1993) 339; 
D.J. Broadhurst and A.L. Kataev, Phys. Lett. 
B315 (1993) 179. 

13. G. Altarelli, P. Nason, and G. Ridolfi, Z. 
Phys. C68 (1995) 257; C.N. Lovett-Turner 
and C.G. Maxwell, Nucl. Phys. B452, 188 
(1995); P. Ball, M. Beneke and V.M. Braun, 
Nucl. Phys. B452 (1995) 563; M. Neubert, 



